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Abstract. Asymptotic properties of random regular graphs are object of 
extensive study in mathematics. In this note we argue, based on theory of spin 
glasses, that in random regular graphs the maximum cut size asymptotically 
equals the number of edges in the graph minus the minimum bisection size. 
Maximum cut and minimal bisection are two famous NP-complete problems with 
no known general relation between them, hence our conjecture is a surprising 
property of random regular graphs. We further support the conjecture with 
numerical simulations. A rigorous proof of this relation is obviously a challenge. 



1. Introduction 

Maximum cut and minimal bisection are two famous problems in graph theory. Given 
a graph, i. e. a set of nodes V and a set of edges E, the goal in the maximum cut 
problem is to split the set of nodes into two groups in such a way that the number 
of edges connecting the two groups is the largest possible. In the minimal bisection 
problem the goal is to split the set of nodes into two equally sized groups in such a way 
that the number of edges between the two groups is the smallest possible. Minimal 
bisection is also known under the name of graph bi-partitioning. Both these problems 
are recognized as NP-complete [l], and both have a large number of appHcations in 
computer science and engineering. Some intensively studied appHcations of the graph 
partitioning problem are circuits design [2j, or data clustering and load balancing in 
parallel computing [3]. For applications of the max-cut problem see, e. g., the survey 
article in Ref. [4]. 

Random r-regular graphs are randomly chosen from all those graphs having N 
nodes and the degree of each node fixed to r. Determining the asymptotic size 
of their max-cut or their min-bisection (bisection width) are classical problems in 
random graph theory, see Refs. O [H [7l HI [9l [TOj for the best known lower and upper 
bounds. However, as far as we know, no explicit relation between the max-cut size and 
bisection width is known in the graph theoretical literature. An exception is provided 
in Ref. [lO], where the same approximative algorithm is used to provide an upper 
bound for bisection width and lower bound for max-cut. 

The main purpose of this note is to conjecture that in random regular graphs the 
size of the max-cut is asymptotically equal to the number of edges minus the size of 
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Figure 1. Two different drawings of tiie same randomly generated 3-regular 
grapti witti A'' = 32 nodes. Left: Example of a minimal bisection of the graph; 
only 6 edges are present between the group of blue (up) and red (down) nodes. 
Right: Example of a maximum cut, only 5 edges are present between two nodes 
of the same color. 



the min-bisection. Our conjecture states that in the large N Hmit, 

\A'IC\ = \E\ - \BW\ + o{\BW\) , (1) 

where is the total number of edges. 

In Fig. [H we present two different drawings of the same 3-regular graph with 
TV = 32 nodes. The left side is a minimal bisection of size \BW\ = 6, the right side 
shows a maximum cut of size \MC\ — rN/2 — 5. It illustrates that Eq. is a highly 
non- intuitive result, since on a given graph there is no straightforward relation between 
the set of edges in the maximum cut and the minimal bisection. As we will show in 
the rest of this note, hints of this conjecture already appeared in various forms in the 
spin glass literature. Our goal is to collect arguments for its justification, state them 
in a language that does not require knowledge of the replica or cavity computations, 
provide evidence from precise numerical simulations, and most importantly, clarify 
the conditions under which this conjecture holds and discuss its generalizations. 



2. Statistical physics formulation of the problems 

In statistical physics, max-cut and bi-partitioning can be formulated in terms of finding 
the ground state of an Ising model on random r-regular graphs. For any graph, the 
general Ising model Hamiltonian reads 

^ = ~ X! Jij^i^j ' (2) 

(i])eE 

where the sum extends over all edges in the graph, Jij is the interaction strength, and 
Si € {—1,-1-1} are the Ising spin variables. The max-cut problem is cast as a ground 
state of the anti-ferromagnetic Ising model, i. e. minimization of ([2]) with Jij = —1 for 
all {ij) e E with respect to the values of the spins {Si}. The min-bisection, or graph 
bi-partitioning, is a ground state of the ferromagnetic Ising model with magnetization 
fixed to zero, i. e. minimization of ^ with Jij — 1 for all {ij) G E subject to a 
constraint ^^Si — 0. Let EQc^{{Jij}) be the energy of the corresponding ground 
state, then the size of the bisection width and the max-cut are 

IBWHM±^2^(iM), |MC| = M_^2MiM), (3) 
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where \E\ is again the total number of edges. 

One can interpolate between the max-cut and min-bisection problems by taking 
the interactions Jij uniformly at random from a bimodal distribution 

P{J.,j) = pS{J., + 1) + (1 - p)SiJ^, - 1) (4) 

and by fixing the magnetization to zero when needed. The disorder in the interactions 
induces frustration on any kind of loopy lattice, in which case the Hamiltonian ^ 
then provides a model for a spin glass [TT]. For sparse random regular graphs the 
conjecture discussed here can be generaHzed as: The ground state energy of ^ is 
asymptotically independent of p, i. e. in the large N Hmit it is 

^Gs({^^,},p) = ^eGs + o(^), {N ~^cx,,0<p<l). (5) 
3. Previous results 

In statistical physics of disordered systems, the replica or cavity method [TTl,fT2] and 
a repHca symmetry breaking scheme ytSj is used to compute the exact ground state 
of Hamiltonian ^ at zero magnetization. Unfortunately, these techniques are not 
rigorous, although in many models their result have been proven, see e. g. Refs. |14[ll5j. 

Using the replica method, Fu and Anderson [16] computed the ground state energy 
of graph bi-partitioning {p = 0) on dense random graphs, i. e. when the degree r = pN 
(0 < p < 1) is a constant independent of the graph size N. Their result reads 

E^T" = UskN'^ Vp{T^) + o{Nh (6) 
where Usk is the ground state energy density of the Sherrington-Kirkpatrick model 
[l7j computed by the Parisi formula [13j, a numerical evaluation giving Usk = 
—0.763219 . . .. They obtained this result by reahzing that on the dense graphs the 
replica equations for Hamiltonian ([2]) at zero magnetization are basically identical to 
the replica equations for the Sherrington-Kirkpatrick model. Moreover, the minimal 
bisection of a graph G plus the maximum-cut of the complement of G (i. e., the graph 
composed of edges that are not present in G) equals (iV/2)^. Using this identity plus 
Eqs. ^ and we obtain that the bisection width is the number of edges minus the 
size of the max-cut |BW| = ~ |MC| + o{N^^^). 

Similarly, using results on the ground state energy of the spin glass, p = 1/2, 
Refs. [I8 t fTO l [20 t I21j computed the bisection width on sparse random regular graphs. 
In the sparse case, the size of the bisection width is linear in the size of the system 
and one therefore obtains |BW| — \E\ — |MC| -I- o{N). This relation can be proven 
rigorously on sparse random graphs with large degree in the first two orders in the 
degree, in particular, \BW\ = rN/2 + UskN^/t -f o{^/r) + o{N) [22] . 

However, the existing literature never discusses for what ensembles of random 
graphs the above results hold. A counter-example is provided by the Erdos-Renyi 
random graphs, where every edge is present with probability a/{N—l). On the Erdos- 
Renyi random graphs the spin glass model, p = 1/2, and the max-cut, p = 1, have 
a positive ground state energy above the percolation threshold, a > 1. In contrast, 
the bisection width, p = 0, of an Erdos-Renyi graph is positive only above a = 2 In 2, 
at which point the giant component reaches size N/2. Thus, we need to discuss 
the theoretical arguments for sparse random graphs and specify the conditions under 
which the max-cut and bisection width are related. 

Note also that on dense graphs, although conjecture ^ holds, the ground state 
energy of ^ is not p- independent. Hence, generalization |(5]) holds only on sparse 
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graphs. E. g., for p = 1 and p = the ground state energy is zero, whereas for p — 1/2 
the ground state energy is UskN^^^- 

For completeness, let us note that the model ^ on random graphs without the 
condition on zero magnetization was studied in Ref . [23] . It was found that there is an 
r-dependent critical value of < pdr) < 1/2 such that for p > Pc{r) the model with 
zero magnetization or non-fixed magnetization are asymptotically equivalent, and for 
p < Pc{r) the two are different as the second develops a non-zero magnetization. From 
the one-step replica symmetry breaking solution, Ref. [23] found, e. g., Pc(3) — 0.142 
and Pc{r) 1/2 for r — > oo. 

4. Theoretical arguments 

We will argue that on sparse random graphs the ground state of ^ at zero 
magnetization does not depend on the fraction of anti-ferromagnetic bonds p. The 
main part of the argument is based on the fact that sparse random graphs (with finite 
mean of the degree distribution) are locally tree like, i. e. the length of the shortest 
cycle passing trough a random node diverges when — > oo. On a tree, all dependence 
on p can be "pushed" to the boundary conditions by using recursively from the root 
the gauge transformations Jij cnJijcrj and Si OiSi, where Ui £ {±1} are chosen 
in such a way as to yield, say, Jij = —1 for all (ij jl. 

Note that this gauge transformation always conserves the Hamiltonian ([2|). So 
we "only" need to discuss the dependence on the boundary conditions. There is a set 
of properties on the boundary conditions that do influence behavior in the bulk of the 
tree, let us call these the relevant properties. To make the connection between the 
system on trees and on the random graph, we need to consider boundary conditions 
on the tree having the same relevant properties as a configuration taken uniformly 
at random from the zero-temperature Boltzmann measure associated with |(2| on the 
random graph. 

Moreover, in order to argue in favor of our conjecture |[5]), we need the 
relevant properties of the boundary conditions to stay unchanged after the gauge 
transformation. A particular relevant property is the total magnetization on the 
boundary conditions. Note that gauge-fiipping of any finite fraction of random bonds 
in the tree causes a random half of the spins on the boundary conditions to flip. Hence, 
only the zero value of magnetization can be treated this way. This is a flrst important 
Hmitation of the conjecture ^ - the p-independence of the ground state of ^ at zero 
magnetization does not generalize to non-zero values. 

We said that all the relevant properties, not only the magnetization, of the 
boundary conditions need to be conserved by the gauge transform. This entails an 
impasse in the mathematical rigor of our discussion and we have to resort to non- 
rigorous arguments impHcit in the cavity method. In the cavity method approach 
Mezard and Parisi [T2| argue that the space of configurations and boundary conditions 
can be split into states. Every state has an associated set of boundary conditions in 
such a way that within each there is no dependence of the bulk properties on the precise 
boundary conditions corresponding to the state. This notion is familiar from the Ising 
ferromagnet in the low temperature phase, where there are two such states. Mezard 
and Parisi [12] treat the case where the number of states grows exponentially with the 
size of the system, this is called replica symmetry breaking. If there is independence 

X Any other pre-defined configuration of Jij's would do as well. 
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of the bulk on the boundary conditions, then there are no relevant properties, and an 
empty set is certainly conserved by the above gauge transformation. Thus, to finish 
our argument we "only" need to show that the solution of the cavity equations (that 
describe the splitting into states) is p-independent or, in other words, conserved by 
the gauge transformation. 

The cavity equations are written in terms of local magnetic fields and their 
distributions over the graph edges (and over the different states, if the corresponding 
problem is glassy j^. It follows from the cavity equations that, if there is a global 
symmetry between positive and negative fields, then the system has to have zero 
magnetization. However, the opposite is not true: requirement of zero magnetization 
does now imply the distribution of fields to be symmetric around zero. The 
inhomogeneity in the graph degree may lead to zero magnetization without overall 
plus-minus symmetry, as can be illustrated again by the example of Erdos-Renyi 
graphs above the percolation threshold, see Ref. [H]. There, at p = the denser 
parts of the graph are more likely to have positive (or negative) fields, and sparser 
parts have excessive negative (or positive) fields. Again, such a non-trivial symmetry 
breaking is not conserved by the gauge transform and, hence, on Erdos-Renyi graphs 
the ground state of ^ may be (and in fact is) p-dependent. Other ensembles of 
non-regular random graphs may also have this degree-fluctuations driven symmetry 
breaking and hence no reason for validity of ([5]) , inparticular for low values of p where 
the equilibrium value of magnetization is not zerqiJ. 

Random regular graphs, on the other hand, have no inhomogeneity in degree 
and they locally look the same from any node in the graph. Moreover, for p = 
and p = 1 there is no inhomogeneity in the interactions Jy either, hence, the cavity 
fields (or their distribution over states) have to be the same on every edge. In such 
a case, the only way to obtain zero magnetization is to have cavity-field distributions 
symmetric around zero. Hence one obtains the same cavity equations for both, the 
graph bisection (p = 0) and the max-cut (p = 1) problems. For the remaining values 
of < p < 1, the neighborhood of every node is different in terms of the set of 
interactions Jy . However, this difference can be pushed to the boundary conditions 
via the gauge transformation. And from the independence on boundary conditions in 
every state it follows that the distribution of fields is the same on every edge. 

In conclusion, the solution of the cavity equations for the ground state of ^ at 
zero magnetization are the same for every < p < 1; this is true on any level of replica 
symmetry breaking and, consequently, the ground state energy of ^ is p-independent, 
as long as the graph of interactions is regular. 

5. Numerical evidence 

We use the extremal optimization (EO) heuristics [251 ES] to find ground states of 
|(2]) at zero magnetization for different values of p. The EO heuristics has been used 
previously for finding ground states on random graphs for graph bi-partitioning (p — 0) 
[Ml [23 [Ml HQ], and for spin glasses (p = 1/2) [;3ll[33]. Thus, it is perfectly suited to 
approximate ground states over the entire range of < p < 1. 

§ In computer science the local magnetic fields are known as the beliefs in the belief propagation 
algorithm, and different states correspond to different belief propagation fixed points. 
II On the other hand the p-independence may be valid for p above some critical value. For example, 
the spin glass p = 1/2 is equivalent to the anti-ferromagnet p = 1 on other ensembles of sparse 
random graphs where the bisection p = is not. 
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A detailed study of the r-EO algorithm in its appHcation to graph bi-partitioning 
and spin glasses is already provided in Refs. |26l [23, and we add only a number 
of minor modifications here. EO considers each vertex of a graph as an individual 
variable with its own fitness parameter. In the graph bi-partitioning, or for any other 

< p < 1, it assigns to each vertex i a "fitness" \i — —bi, where bi is the number 
of "bad" (unsatisfied) edges connecting i to other vertices. At all times an ordered 
list is maintained, in the form of a permutation 11 of the vertex labels i, such that 
An(i) < An(2) < • ■ ■ < An(Ar), and i = n(fc) is the label of the k-th ranked vertex in the 
list. In its most elementary version, EO forces sequential updates of the momentary 
worst variable i — 11(1) at any update step, irrespective of the outcome, inducing 
a cascade of adaptive reorderings in the list. Since all variables occupy an identical 
and 0(l)-sized state space, Aj = 0,-1, . . . , — r, for r-regular graphs, the list is highly 
degenerate and maintaining order or selecting variables (with fair tie-breaking rules) 
is done in 0(1) computations. 

To define a local search of the configuration space, we must define a 
"neighborhood" for each configuration within this space. At zero magnetization for 
all < P < 1, as an improvement over our previous implementation of EO for 
graph bi-partitioning, we proceed here by allowing imbalanced partitions up to a 
margin of ±2 vertices, independent of system size N. Then, we can pursue single- 
variable updates as long as the resulting configuration remains within the allowed 
imbalance. Valid ground states are only accepted, if the current partition is perfectly 
balanced (although any 0(1) imbalance for increasing N should result in identical 
scaling behavior). In this form, the single-flip neighborhood trivially generalizes to 
spin glasses at freely fluctuating magnetization, for which we simply ignore whether 
partitions remain balanced within the margins or not. 

Much improved results are obtained with the following one-parameter 
implementation [25], called t-EO: An integer 1 < fc < iV is drawn from a probability 
distribution P{k) oc fc~'^, 1 < k < N, on each update, for flxed r. Then, the vertex 

1 = n(fc) from the rank-ordered list of fltnesses is selected for an unconditional update. 
Over the course of a run (here with tmax = O.liV'^ update steps), the costs of the 
configurations explored varies widely, since each update can result in better or worse 
fitnesses. The cost minimum of the best configuration seen during all runs for an 
instance is the output of the EO-algorithm. We choose at least three uncorrelated 
restarts for each instance here. The number of restarts is automatically adjusted for 
each instance such that twice as many runs are undertaken than was necessary to 
encounter the putative ground state for the first time. Only at larger system sizes and 
degree r, when more than 10% of instances require more than those three restarts, we 
initially set the duration of each run to up to tmax = 0-5N^ update steps. 

Note that no scales to limit fluctuations are introduced into the process, since the 
selection follows the scale-free power-law distribution over ranks P{k) and since all 
moves are accepted. Instead of a global cost function, the rank-ordered list of fitnesses 
provides the information about optimal configurations. This information emerges in a 
self-organized manner, merely by selecting with a bias against badly adapted variables, 
rather than ever "breeding" better ones. A theoretical analysis of the optimal r- value 
is discussed at length in Refs. [Ml EH IS] , here, some initial trials suggest optimal 
values of r = 1.2 — 1.3, which we have used throughout. 

Let us call ecs {p, N) the ground state energy density, averaged over graphs and 
disorder in interactions, of ^ at magnetization fixed to zero, with p being the fraction 
of anti-ferromagnetic edges and N the graph size. Denote by eGs{p,N) the same 



Conjecture on the maximum cut and bisection width in random regular graphs 



7 



quantity at an arbitrary magnetization. We have obtained ecsip, N) on random 
regular graphs of degrees r between 3 and 10, and graph sizes between TV = 32 and 
1024. Statistical errors of our averages have been kept small by generating a large 
number of instances for each N and r, typically n/ ~ 10^ for N < 200, nj « 10^ for 
N > 256. 

All our data are indeed consistent with the conjecture that on sparse random 
regular graphs the ground state energy of ^ at zero magnetization is p-independent 
in the thermodynamic limit, N ^ oo. However, we also observe that the finite-size 
correction are p-dependent. This can be understood intuitively on the very particular 
case of random 2-regular graphs. A random two-regular graph is basically a set of 
cycles of length ~ logiV. Hence, whereas the bisection width is at most 2 edges, the 
number of edges minus the max-cut size is of order N/ log N (a unit cost for every 
other cycle in the graph). For r > 3, the finite size correction are not that large, but 
they are quite different for different values of p. 



r 


eiRSB 


ecs [p = 5) 


ecs {p = 0) 


c^(p = 0) 


3 


-1.27231 


-1.2716(1) 


-1.2704(2) 


0.89 


4 


-1.47295 


-1.472(1) 


-1.469(1) 


0.92 


5 


-1.67520 


-1.673(1) 


-1.6717(5) 


0.86 


6 


-1.82917 


-1.826(1) 


-1.824(1) 


0.87 


7 


-1.99566 


-1.990(3) 


-1.990(1) 


0.85 


8 


-2.12681 


-2.121(1) 


-2.120(2) 


0.87 


9 


-2.27093 


-2.2645(5) 


-2.263(3) 


0.85 


10 


-2.38769 


-2.378(3) 


-2.379(4) 


0.86 



Table 1. Asymptotic ground state energy per spin for different values of the 
graph degree r. The second column, euigB, present p-independent one-step 
replica symmetry breaking results [35 j . The third column, eGs(l/2), contains the 
numerical values of the extrapolated ground state energy for the spin glass [3 1|. 
p = 1/2. The fourth and fifth column give ground state energies eGs(O) and 
scaling coefficients i^{0) for the graph bisection problem obtained from infinite 
graph-size extrapolations according to Ul, as shown in Fig. [21 With minor 
exceptions for the smallest degrees, ecs (1/2) and eGs(O) are the same within 
error bars. Moreover, lacking a theoretical justification for Eq. JTJ, the effective 
error bars of the fitted values are larger than denoted in the table. 

In Tab.[ll we compare the asymptotic ground state energy densities for different 
values of graph degree r. The second column in this table presents the p-independent 
one-step replica symmetry breaking results for the ground state energy of ([2]) at zero 
magnetization, computed with the formalism developed in Ref. [35]. Note that the 
exact value for the ground state would be provided by the full-step replica symmetry 
breaking scheme which would give slightly larger values. The data in the third column 
are taken from Ref. [31]. In this case the magnetization was not fixed to be exactly 
zero on every instance, but it is zero in density in the thermodynamic limit. The data 
in [3T| [32] are consistent with a power law scaling 

eosiN) = eos + aN-^ . (7) 

with the value of the exponent a; = 2/3 for all r. The finite-size scaling for the graph 
bisection is clearly not consistent with uj = 2/3, as is illustrated in Fig. [2] separately 
for odd and even values of As was noted in Refs. [311 133) . ground state energies 

1 Note, however, that due to possibly strong higher order corrections to jTj the values of a; in Table[T] 
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are strongly affected by the existence or absence of "free spins" on graphs with purely 
even or odd degrees, respectively. 
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Figure 2. Plot of the rescaled average ground-state energy densities, eQs{l^)/V^, 
for the bi-partitioning of random regular graph of degree r as a function of . 
Errors are smaller than symbol sizes and are omitted for clarity. Non-linear fits 
to each data set according to Eq. JtI are indicated by dashed lines. Symbols on 
the ordinate (green diamond for odd r, red star for even r) mark the extrapolated 
values from Table [T] (rescaled by l/\/r) of the corresponding spin glass ground 
states from Ref. [31]. For large even and odd degree r, the extrapolated values 
approach the ground state energy of the Sherrington-Kirkpatrick model UsK 
(black cross). The data are consistent with the conjecture that the spin glass 
asymptotic ground states are equal to the graph bi-partitioning ones. 

On the left hand side of Fig. O we study numerically the dependence of average 
ground-state energies as a function of p at finite sizes on 3-regular graphs only. While 
there are significant differences in the magnitude of corrections - even on average - 
for the smaller sizes, finite-size behavior soon becomes virtually independent of p and 
appears to converge towards the thermodynamic values found consistently at p = 
and 1/2 listed in Tab.[Tl 

On the right hand side of Fig. [Sj we address the question on how finite-size 
corrections are impacted by the constraint on magnetization being fixed to zero. 
We compare the average ground state energy of the anti-ferromagnet (p = 1) with 
magnetization strictly zero and with no constraint on magnetization. Clearly, in the 
unconstrained case the distribution of ground-state magnetizations is symmetrical 
above pc, with vanishing fiuctuations in the thermodynamic limit. Hence, the 
average magnetization must be zero. Yet, at finite size, it seems conceivable that 
fixing the magnetization makes otherwise lower-energy states, i. e. ground states 
of the corresponding unconstraint model, unattainable, which may shift the ground 
state energies of the constraint system upwards. Indeed, as Fig. [3] demonstrates, 
while indistinguishable thermodynamically, average energies increase by large amounts 
relative to the unconstraint case especially for small sizes. But the actual finite-size 
corrections seems to be compatible with uj = 2/3 in both cases. In our simulations, 
we generate r-regular graphs in such a way that multiple edges between identical 
vertices are not forbidden in their random assignment. Such multi- linked vertices 
have a probability of ~ 1/N and, although noticeable at small size, do not affect any 
asymptotic scaling, see Fig. [H Forbidding such edges makes it difficult to generate 
valid graphs especially at larger r and small sizes. 



may be skewed. 



Conjecture on the maximum cut and bisection width in random regular graphs 



9 



r = 3 



•- -• N= 32 

-- N= 54 
•- -■ N= SO 

N=12S 
N=160 
» -« N=256 



-1.18 
-1.20 
Z -1.22- 

-1.26 - 
-1.28 - 



r = 3 (p=l) 



0.2 0.3 0.4 0.5 0.6 
P 



0.7 0.8 0.9 







• M = 

o M unconstraint 

o M uncon. (True Bethe) 



1/N' 



0.06 
2/3 



0.08 



IT 



0.1 



Figure 3. Left: The average ground state energy density for finite-size 3-regular 
graplis as a function of the density of anti-ferromagnetic bonds, p. The red dashed 
line marks the value extrapolated for p = or 1/2 from Tab. [T] The data are 
consistent with the conjecture that the ground state energy is asymptotically 
independent of p. Right: Plot of average ground state energies for the purely 
anti-ferromagnetic spin model (p = 1) on 3-regular graphs as a function of system 
size, both, at fixed (M = 0) and at unconstraint magnetization. The constraint 
case appears to extrapolate well linearly on a N~^^^ scale without any transient, 
while the unconstraint case shows a small deviation of about 5% when all data 
is fitted according to Eq. ([7|- With blue squares, we also plot data for the same 
problem with unconstraint M but disallowing multiple edges between two vertices. 
All other data was obtained allowing graphs with such edges. 



6. Possible generalizations 

As we have argued above, the relation between max-cut and bisection width does not 
generaUze (at least not straightforwardly) to the case of non-zero magnetization (not 
equally-sized groups) nor to the case of non- regular graphs. However, it does generalize 
to finite temperature properties of the Hamiltonian ^ at zero magnetization. 

The relation between max-cut being the number of edges minus the bisection 
width also generalizes to the case of hyper-graphs. In statistical physics, corresponding 
models are know as models with p-spin interactions, in computer science as the XOR- 
SAT problem (boolean constraint satisfaction problem consisting of sets of linear 
equations) . 

An alluring but not (fully) valid generalization to discuss is the case of Potts spins 
Si = 1, . . . , fc. The max-cut problem is then replaced by max-fc-coloring problem and 
the bisection by fc-partitioning. According to a result by Kanter and Sompolinsky 
[36], analogous to the one of [16], in dense graphs the two problems are related. Be 
pN the degree of the graph, and k the number of colors, then according to Ref. [36] 
the maximum number of non-monochromatic edges is 

iMaxColl = ^N'[l - i) + TvlM^RT^, (8) 

whereas the minimal number of edges between groups in the best balanced fc-partition 
is 

|fc-part| = |7V2(l - i) - iviM^^(r3^. (9) 

Here, U (k) is in both the expressions for the ground state energy of the fully-connected 
Potts model with k colors [37], numerically given in Ref. [36], while in the large-fc limit 
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it is limfe^oo U (fc) = y/ (kink). For fc = 2 the relations (UHll) reduce to the Fu and 
Anderson result 

Based on the analogy of the dense graph case, we would thus expect a 
generalization for sparse graphs, p = c/N, also for k > 2. However, for three and 
more colors there is no apparent relation between the max-coloring and ^-partitioning. 
Whereas there are some version of the Potts glass equivalent to the coloring problem, 
see e. g. Ref. [38], there is no obvious Gauge transform able to transfer the Potts 
ferromagnet (allowing 1 out of q values) on the Potts anti-ferromagnet (allowing q—1 
out of q values). This can be seen explicitly in the difference between the replica 
symmetric equations for the two problems. In the warning-propagation sense [39] . 
the neutral warning in max-coloring is created, if two colors have the same value of 
an incoming field and the third one has a larger value. In contrast, in partitioning 
the third value needs to be smaller. Also, the expressions for the repHca symmetric 
energy are different, even after numerical evaluation. The equivalence thus holds only 
asymptotically in the first two orders of the degree of the graph (as suggested by the 
result in Ref. [Mj). This underlines the exceptional nature of our main conjecture ([T]) 
for /c = 2. 

7. Conclusion 

In this note we describe, explain, and support by numerical evidence a conjecture 
that on sparse random regular graphs the ground state energy value of the spin glass 
Hamiltonian ^ at magnetization fixed to zero does not depend on the fraction of 
anti-ferromagnetic bonds. Although hints towards this conjecture can be found in the 
existing literature, we state it as a clear mathematical conjecture understandable for 
non-specialist in spin glass theory: In random regular graphs, the asymptotic size of 
the max-cut equals the number of edges minus the minimal bisection width. We also 
summarize necessary conditions and limitations of this conjecture, in particular, that 
it does not generalize (at least not in a way we could see) to non-regular graphs and 
non-zero values of the magnetization. We also support the conjecture by extensive 
numerical evaluations of the ground states. Finally, we are positive that this note will 
be useful for the mathematical and computer science community and that it will lead 
to a proof of this conjecture in the near future. 
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